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We discuss possibility of the existence of tetraquark states made of four quarks in the string (flux 
tube) model. The new particle is composed of a diquark and an anti-diquark which are connected by 
a color flux. It is shown that the vibrational and rotational excited states of the string explain some 
non-qq mesons observed experimentally. Moreover we discuss the decay widths of such tetraquarks 
with the use of the Schwinger mechanism. 
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I. INTRODUCTION 

It is well known that there are only two kinds of 
hadrons in Nature: mesons and baryons. The former are 
made of a quark and an antiquark and the latter made of 
three quarks. Their existence is explained by quantum 
chromodynamics (QCD) 1] that is the fundamental the- 
ory for hadron physics. There is, however, a possibility of 
the existence of other type hadrons, such as four-quark 
particles. In fact a possible four-quark state has been 
suggested experimentally by the several groups Q. 

In this paper we discuss the possible existence of a 
four-quark state which we call a "tetraquark" ; it is made 
up of two quarks and two antiquarks. Since Jaffe pre- 
dicted such hadrons with the use of the MIT bag model 
0, 3 , many non-qq mesons have been discussed by many 
authors Particularly it is conjectured that the 

nonet of light scalar mesons [/ (600), / (980), a (980), 
k(800?)] with masses below 1 GeV seem to fit a qqqq 
description. Instead a heavier scalar nonet above 1 GeV 
seems to be largely qq configuration. The light scalars 
seem to have a significant component of qqqq and the 
pairs of quarks (and antiquarks) are separately corre- 
lated in a flavor, color and spin antisymmetric state. A 
diquark and an anti-diquark neatly couple together to 
give a nonet of color singlet scalar mesons. There is also 
evidence for strong diquark correlations in the baryon 
spectrum 

We will consider a possibility that this strong di- 
quark correlation leads to a molecular-like configuration, 
(qq) — (qq) with the use of the chemical bond-like ex- 
pression. To this end, we will study the four-quark sys- 
tem by using the hadron string model or the flux tube 
model. According to the model, mesons are made up of 
a quark and an antiquark which are connected by the 
classical gluon field 13]. The baryons are also made up 
of a quark and a diquark connected by the same flux 
Using this simple model, F.Takagi and one of 
the present authors (M.I.) have shown that many excited 
states of hadron observed experimentally can be identi- 
fied as the vibrational and rotational excited states of 
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this one-dimensional string 0, . 

Since any hadron must be a color-singlet, the only pos- 
sible four-quark state is a qqqq state. There are here two 
configurations: (a) (qq) — (qq) and (b) (qq) — (qq)- We sup- 
pose that other configurations (q) — (qqq) and (q) — (qqq) 
are not so stable. Assuming the one-gluon exchange in- 
teraction, the qq system is attractive in the color singlet 
state (meson) so that the system (a) probably decays into 
two mesons at once. On the other hand qq (qq) system 
is attractive in the color anti-triplet (3bar) state, which 
means that it is equivalent to q (q) in the color SU(3) 
space. That is the system (b) is the same as the me- 
son as far as the color degree of freedom. Therefore a 
four-quark state considered in this paper is assumed to 
be composed of a diquark and an anti-diquark connected 
by the same gluon flux as the meson. 

It should be emphasized here that our approach is dif- 
ferent from that in Jaffe's papers. According to them, 
the nonet of light scalar is described by the s-wave qqqq 
sector in the MIT bag model; all the quarks are in the 
same s-state in the bag potential. On the other hand, a 
scalar in our model is described by the s-sate of the di- 
quark and anti-diquark where each quark is also s-state. 
The latter s-state is quite different from the former one. 
The wave function of the quark would be a superposi- 
tion of many states of the higher angular momenta in 
the MIT bag potential. The relation of our model to the 
MIT model is analogous to that of the cluster model to 
the shell model in nuclear structure models. Therefore 
the tetraquark discussed in this paper are not the nonet 
of light scalar but other heavier mesons. 

In the next section (Sec. II), we present our flux tube 
model and derive the mass spectra and decay widths us- 
ing the WKB approximation and the Schwinger mecha- 
nism. The numerical calculations of these quantities are 
carried out in the Sec. III. Finally (Sec. IV) several discus- 
sions and summary are given. 



II. A FLUX TUBE MODEL 

According to [l^ , our system is described by the rela- 
tive distance r between the diquark and the anti-diquark 
and the rotation angle 9 of the string. The geometry is 
shown in Fig.l where their masses is denoted by Mi and 
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FIG. 1: The schematic picture of the hadron string. 



Mi- We take the origin as the center of mass. Since 
our discussion is restricted to the nonflavored mesons, 
we assume M\ = M2 = M. Then the Lagrangian of this 
system is given by 



L{r, r, w) = -2M \/l - v 2 - V(r, u), 



(1) 



where v = (r +r u> )/4 is the square velocity of each di- 
quark. The first term represents the action of the diquark 
and the anti-diquark and the second one the potential en- 
ergy of the system which is represented by 



V(r,u>) = 2 



V2 



iy 1 — r 2 u> 2 dr 



(2) 



The first term denotes the energy of the gluon flux with 
the string tension a (the energy per unit length) [Isf 
and Vl — r 2 uj 2 comes from the Lorentz contraction of 
the transversal motion. The second one represents the 
one-gluon exchange potential with the coupling constant 
e. This expression for the potential is known as Cornel 
potential and confirmed by the Lattice QCD simulation 
[T9L l20| . It should be noted that this Lagrangian is also 
valid for the qq system (ordinary mesons) because their 
color charges are the same as those of the qq — qq one. 

The Hamiltonian of our system is defined by H = rp + 
u>l — L where p = dL/dr and I — dL/d6. Then it is 
written as 



H(p,l,r) = 



2M 



r/2 



2ft 



Vl — T 2 U1 2 



zdr . 



(3) 



where the factor Vl — r 2 u> 2 in the right-hand side repre- 
sents the Lorentz contraction along the transverse direc- 
tion of the string. The excited states of the system are de- 
scribed by the Schrodinger equation, H(p, I, r)ip(r, 9) — 
Ei/j(r,8). Carrying on the separation of variables tp = 
(x/r) expimO, this equation is reduced to 



U (' -i-^-,m,r 



X(r) = Ex(r), 



(4) 



where I is the angular momentum quantized as m — 
0, 1, 2, 3, • • •. This equation can be solved using the WKB 
approximation and the excited energy E is obtained by 
the Bohr-Sommcrfcld formula, 



p r dr = (n + l/2)ft. (n = 0, 1, 2, 3, • • •) (5) 



The right hand side should be replaced by (n + 3/4) h 
when to = 0. The full derivation of the formula is ex- 
plained in Ref. 

Here let us discuss the qualitative properties of the 
mass spectra (neglect the masses and the one-gluon ex- 
change term for convenience). In the case of to = (vi- 
brational motion), the energy becomes E = 2p+ar = aro 
(r is the maximum length). This leads to <f pdr — 
(1/2) arg 2-7m. Thus we have 



E = Aiian. 



(6) 



In the case of n — (rotational motion), the energy is 
E = 2m /r + ar if the kinetic energy is approximated 
b y the c entrifugal potential. This is minimized by r — 
\J 2tti I a so that E = 2ar. These relations lead to 



E z = 8am. 



(7) 



In both cases, it is concluded that the energy squared 
E 2 is proportional to the quantum numbers n or m and 
the ratio of the proportional constants is approximately 
7T : 2. 

Next let us consider the decay width of our system. 
The possible decay process that we consider is a pair 
production of qq inside the flux tube which is known as 
Schwinger mechanism j2^|-[2^|. The probability of the 
pair production in a unit space-time volume in the tube 
is given by 



w 



2 00 1 
a \ - x - 1 

4-7J-3 Z-^/ Z-^/ r,2 



exp 



q n—1 



(8) 



where a is the string tension and indicates a sum- 
mation over all quark flavors with mass m q (q=u,d,s) 
[23|- H3- The probability of the decay (pair production) 
in the time interval dt is given by 



dW 



\(r) | 2 wv(r)drdt, 



(9) 



where v(r) denotes the volume of the string with the 
length r and is given by 



r/2 



v(r) = 2a \J 1 — r 2 uj 2 dr. 
Jo 



(10) 



In the right-hand side, a denotes the area of the cross 
section of the tube in the rotating coordinates system. 
The decay width of our system is given by 



T = = I \x(r)\ 2 WVlr)dr. 



(11) 
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FIG. 2: The vibrational excited states of the rho meson. The 
quantum number of the vibrational motion is denoted by n. 



FIG. 3: The vibrational excited states of the tetraquark. The 
quantum number of the vibrational motion is denoted by n. 



This formula gives the decay widths of the excited 
hadrons due to the quark pair production. As is shown 
in this formula, the decay width is proportional to the 
length of the string. This means that it increases as the 
excited energy increases. 



III. NUMERICAL RESULTS 

We are in a position to carry out the numerical calcu- 
lations of our tetraquark particles. Our string model has 
two components, the color flux and the diquark. The for- 
mer is characterized by the string tension a and the cross 
sectional area a and the latter is concerned to the diquark 
mass. The coupling constant e will be used to be e = 0.4, 
which is a usual value 0] . The parameters concerning to 
the color flux should be taken as the same values as those 
of q — q system, because a diquark (anti-diquark) has the 
same color charge as that of an antiquark (quark) . These 
parameters have been determined so as to reproduce the 
experimental mass and decay width of p(770). 

When we take a = O.llGeV 2 and (Mi = M 2 = 
lOMeV), the calculated vabrational spectrum of the rho 
meson is shown in Fig. 2. In the case of m = 0, the 
interval of the integral (5) is < r < ro so that the inte- 
gral diverges. Therefore the cutoff parameter e = O.Olfm 
(e < r < r ) has been introduced so as to reproduce the 
mass of p(770). Its decay width (r = 150MeV) deter- 
mines the other parameter, a = 75GeV~ 2 . We use these 
parameters in later calculations hereafter. 

As for the diquark mass, we take the most simple di- 
quark which is composed of u and d quarks. In other 
words we restrict ourselves to the flavor SU(2). In this 
case, the most stable diquark is color anti-triplet, flavor- 
singlet and spin-singlet taking into account the Fermi 
statistics. Since the nucleon is made up with this di- 
quark and the remaining quark, our diquark mass is de- 
termined by the nucleon mass N(939) as is done in Ref. 
0: M — 300MeV. 

First it is adequate to discuss the possible quantum 



numbers of our tetraquarks, which has configuration 
(qq) — (qq). Since the diquark (anti-diquark) is flavor- 
singlet and spin-singlet, it is isoscalar (/ = 0) and has not 
spin angular momentum. The total angular momentum 
J is nothing but the orbital one: J = m. The G-parity 
operation defined by G — Ce l7rl2 is the same as that of 
the ordinary parity P, because the charge conjugation C 
is equivalent to P in our system. Therefore we should 
seek mesons with I G (J PC ) = ± (J ±± ). 

Let us discuss the vibrational excited states with the 
angular momentum m = 0. The square of mass is drawn 
as a function of the quantum number n in the Fig. 3. The 
closed circles in the figure represent the corresponding 
particles with 0+(0++) tentatively. They are /(1110), 
/o(1500), /o(2020), / (2200) 21]. The head state /(1110) 
is + (even ++ ) and its spin J is not determined except for 
J =even integer. It is seen our simple model reproduces 
such particles fairly well. The trajectory is an almost 
straight line as is shown in the Eq.(6). 

Next let us consider the rotational excited states with 
n = 0, which is called Regge trajectory. In the Fig. 4, 
the square of mass is shown as a function of the angular 
momentum. As discussed above, the quantum numbers 
of the Regge trajectory should be + (0 ++ ), _ (1 __ ), 
+ (2 ++ ), CP (3 ), • • •. So the closed circles represent the 
corresponding particles tentatively: /(1110), wi(1420), 
/ 2 (1640), w 3 (1945), / 4 (2050), lu 5 7, / 6 (2510). The candi- 
date for cc>5 has not been observed yet. We see that our 
model explains such particles fairly well. The deviation 
from the straight line is due to the diquark mass. 

Next we go over to the discussion on the decay width 
of the tetraquark. Here it is interesting to notice that our 
tetraquarks with E < 2GeV can not decay with the pro- 
cess of the pair production by virtue of the conservation 
of energy. Because the total energy of the decay products 
(a baryon and an antibaryon) must be larger than about 
2GeV. Therefore it is predicted that the decay widths of 
/(1110), /o(1500) are suppressed and those of / (2020) 
and /o(2200) are enhanced. The results are plotted as 
circles on the (M,T) plane [2l| as shown in Fig. 5. The 
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FIG. 4: The rotational excited states of the tetraquark. The 
orbital angular momentum is denoted by m 



FIG. 6: The vibrational excited states of the q — q mesons. 
The quantum number of the vibrational motion is denoted by 
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FIG. 5: The masses and decay widths of all the unflavored 
mesons observed experimentally are plotted by the open cir- 
cles. The filled ones represent the candidates for the vibra- 
tional excited states of the tetraquark. The corresponding 
calculated states are denoted by the triangles. 



decay mode to baryon-antibaryon has not been observed 
yet in the case of / O (2020) [2^, H3|- The experiments in 
Fig. 5 is focused on the process pp — > pf{Air)p s . We con- 
sider that even if a baryon-antibaryon pair is created in 
the process, the antibaryon would interact at once the 
other baryon and decay into pions. As for the decay of 
/o(2200), it is seen in the experiment pp — > tttt where 
only the 7T7t mode is selected j^. Although it is not 
detected in the experiments, we conjecture a resonance 
reaction of pp — > pp. This resonance would support the 
existence of the tetraquark considered in our model. 

The calculated decay widths are shown by the tri- 
angles. The black circles represent /(1110), /o(1500), 
/o(2020), /o(2200) and the open ones all the other un- 
flavored mesons observed experimentally. The enhance- 
ment of the decay widths of / (2020), /o(2200) may 
be explained by the open of the decay channel due to 
the pair production. The decay widths of /(1110) and 
/o(150Q) are relatively small and caused by the other de- 
cay modes. Of course it is not adequate to expect quan- 



titative discussions on the experimental values by our 
model. Because our model is oversimplified and the data 
on the decay widths have large ambiguity. Here we are 
satisfied with that the calculated decay widths are the 
same order of magnitude as the experimental data. 



IV. DISCUSSIONS AND SUMMARY 

In the previous section, we have studied the scalar 
mesons above lGeV from the standpoint of the qq — qq 
model. There are, however, the other scalar mesons not 
described by the present model. In fact they seem to be 
/o(1370) and /o(1710). There is a possibility that they 
are described by the usual q — q model. We will investi- 
gate such scalar mesons using the simple q — q model in 
this section. 

Let us consider a meson composed of q and q con- 
nected by the color flux. If the meson is p-wave, we have 
the scalar meson, because the intrinsic parity of q is mi- 
nus. The mass spectrum of this system is calculated in 
the same manner as was done in the previous section. 
In particular we restrict ourselves to the case of m = 1 
(p — wave). The vibrational excited states are shown as a 
function of n in Fig. 6. The closed circles in the figure rep- 
resent the corresponding particles /o(1370) and /o(1710) 
tentatively. Therefore it is suggested that these particles 
are composed of a quark- antiquark in the p-wave state. 

What is a discrimination between the qq — qq and q — q 
configurations? One possible clue is the decay width. 
That of the four-quark particle is suppressed for the con- 
servation of energy as mentioned before. In fact the decay 
width of /o(1370) is larger than 1 GeV. However that of 
/o(1710) is not so large, that is 140 MeV, which is larger 
than those of /(1110) and /o(1500) [2l|. 

In conclusion, we have studied possibility of the exis- 
tence of the tetraquark states which are composed of a di- 
quark and an anti-diquark connected by a color flux. The 
vibrational and rotational excited states of the string are 
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calculated with the use of the WKB approximation. It 
is suggested that they may be identified as the observed 
mesons. Moreover the decay widths of the tetraquarks 
are calculated using the Schwinger mechanism. They are 
suppressed by the conservation of energy, which is consis- 
tent with the experimental results. It is a very interesting 
problem to study the possibility of the tetraquark. 
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